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A model with hyperchaos is studied by means of Lyapunov two-parameter analysis. The regions 
of chaos and hyperchaos, as well as autonomous quasiperiodicity are identified. We discuss the 
picture of domains of different regimes in the parameter plane of coupled systems, corresponding to 
the cases of interaction of quasiperiodic and hyperchaotic subsystems. 
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A hyperchaotic attractor is characterized by the pres¬ 
ence of several positive Lyapunov exponents. Dynamics 
of systems with hyperchaos is more complex and diverse 
than that of the system with simple chaos. Examples of 
hyperchaos were pointed out yet by Rossler and others 
researchers [1-6]. Recently some new suitable examples of 
models with hyperchaos were discovered. A study of such 
problems as scenarios of the emergence of hyper chaos, 
structure of the hyperchaotic attractor, a possibility of 
hidden attractors et al. was undertaken in refs. [7-11]. 
It should be noted that such systems are promising from 
the point of view of secure communication because they 
can realize high-dimensional chaos. However, many ques¬ 
tions in the theory of hyperchaos remain open. For ex¬ 
ample, the problem of the interaction of the hyperchaotic 
systems is important and promising. In this report we 
discuss some aspects of this problem. 

Let us consider a model suggested recently in paper 
[7]. It is a modified form of well-known Lorenz system: 

x = a(y — x) + eyz, 
ij = cx - dxz + y + u, , * 

z = xy — bz, ' ' 

u = —ky. 

Here x,y,z,u are the dynamical variables and 
a, 6 , c, d , e, k are control parameters. 

The model © has several fundamental properties and 
advantages [7]: 

1 . it is dissipative at (a + b — 1 ) > 0; 

2. equations © are symmetric with respect to z-axis: 
(x, y, z, u ) ->■ (-X, -y, z, -u); 

3. it has only the trivial zero equilibrium. 

For fixed parameter values a = 35, c = 25. d — 5, e = 
35, k = 100 Chen et al. [7] undertook the one-parameter 
analysis of model (pQ). In Fig.l we reproduced the spec¬ 
trum of Lyapunov exponents and typical phase portraits 
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at varying parameter b. Note that, as usual, one expo¬ 
nent is always zero A 3 = 0 and is not shown in the figure. 
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FIG. 1: Plot of the Lyapunov exponents as a function of 
parameter b for the model (jTj) . (Exponent A3 = 0 is not 
shown in the plot.) 


The system © exhibits hyperchaos in a wide range of 
parameter values. It is interesting that along with such a 
regime, as well as periodic and chaotic regimes, the model 
also demonstrates quasiperiodic dynamics. Thus, model 
© is one more example of an autonomous system with 
quasiperiodicity (see. also [13,14]); this fact increases its 
importance. 

One-parameter analysis does not give us full informa¬ 
tion about typical regimes and their localization. That 
is why we turn to the two-parameter analysis. Let us use 
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charts of Lyapunov exponents [15-17]. We calculate the 
values of the Lyapunov exponents at each point of the 
plane (fc, b ) and color this plane in accordance with the 
type of observed regime: 

• periodic regime P (all exponents are negative); 

• quasiperiodic regime X 2 (one zero exponent); 

• chaos C (one positive exponent); 

• hyperchaos HC (two positive exponents). 

(One trivial zero exponent is not taken into considera¬ 
tion.) The corresponding chart of Lyapunov exponents 
is shown in Figure 2. The color palette is presented to 
the right of the figure. It can be seen that the area of hy¬ 
perchaos entirely occupies extensive part of the plane of 
the control parameters. The domain of existence of two- 
frequency tori is visualized very well. We can distinguish 
the boundary of quasiperiodicity region - the Neimark- 
Sacker bifurcation line with the Arnold tongues emerging 
from it. 



FIG. 2: Chart of Lyapunov exponents of the system (pQ). The 
color corresponds to the type of regime. Two lines are shown 
for which the scan of parameter plane is executed. 


Let us study now dynamics of two coupled models m, 
the coupling terms are represented by the differences of 
the corresponding variables: 

£1 = a(yi - xi) + eyizi + y, (x 2 - xi), 

Vi = cx 1 - dxizi + yi + «i + y (y 2 ~ Vi ), 
h = xiVi - Mi + A* (Z2 - Zi ), 
ui = -hyi + y(u 2 - mi) , 

x 2 = a(y 2 ~ x 2 ) + ey 2 z 2 I // (j't - x 2 ), 
y 2 = cx 2 - dx 2 z 2 + y 2 + u 2 + y (yi - y 2 ) , 

£2 = x 2 y 2 - b 2 z 2 + y(zi~ z 2 ), 
u 2 = -k 2 y 2 + y (ui - u 2 ). 

( 2 ) 


Here x±, y±, Z\U\ and £ 2 , 2 / 2 , £ 2,^2 are dynamical vari¬ 
ables of the first and the second subsystems, respectively, 
and 11 is the coupling parameter. We fixed a set of pa¬ 
rameters a = 35, c = 25, d = 5, e = 35 and change 
coefficients &i, & 2 , 5i, 62 - 

The presence of set of Arnold tongues of different peri¬ 
ods in the autonomous system allows us to associate with 
the model the so-called ”frequency parameter”, corre¬ 
sponding to the motion along the Neimark-Sacker bifur¬ 
cation line. We choose two routes corresponding to vary¬ 
ing of this effective parameter value - one lies through the 
region of quasiperiodic regimes characterized by built-in 
set of Arnold tongues (Route 1 in Figure 2 ) and the sec¬ 
ond through the area of hyperchaos (Route 2 ). 

These routes correspond to variation of the frequency 
parameter in the second subsystem. In the first subsys¬ 
tem let us fix the values k± , b\ so that we can observe ei¬ 
ther quasiperiodicity or hyperchaos. The values of these 
fixed parameters are marked by points on lines in Figure 
2 . 

The first scenario. Let us fix parameters k\ = 430, 
b\ = 12 in the first subsystem and plot a chart of Lya¬ 
punov exponents. Herewith along the y-axis we will vary 
coupling parameter /i and along x-axis - coefficient &2 
measured along the marked Route 1 in Figure 2 . This 
case corresponds to interaction of quasiperiodic oscilla¬ 
tions. The resulting chart of Lyapunov exponents is dis¬ 
played in Figure 3. Now we have an additional possi¬ 
bility of the presence of a three-frequency tori T 3 (two 
zero exponents) and four-frequency tori T 4 (three zero 
exponents). 

At a low coupling parameter of the model (|2j) it is 
naturally to expect an emergence of four-frequency torus 
- this hypothesis is confirmed by the numerical study. 
Also on the chart we can distinguish wide area of two- 
frequency tori in form of peaked tongue touching the line 
of zero coupling. The point of contact corresponds to 
the condition k<i = Aq 5 62 61 when the subsystems 

are identical and their ” frequency parameters” coincide. 
Also, there are narrow tongues of X 3 torus embedded 
in the region of T 4 tori. In right part of Figure 3 we 
demonstrate the phase portraits and Poincare sections at 
points of presence of X 2 , X 3 and T 4 tori. A very narrow 
bands of periodic regimes are visualized. 

Note that inside the area of two-frequency tori it is 
possible to observe a phase synchronization - some kind 
of synchronous quasiperiodicity. Illustrations of such a 
regime are given in Fig.4. We display the phase por¬ 
traits of both partial systems of model ([ 2 |), for which it is 
possible to determine the phase of the oscillations (Fig¬ 
ure 4a,b). Fig.4c represents the time-dependence of the 
phase difference ipi — (f 2 ] phase is calculated according 
to the formula cpi = arctan (yi/xi) ,<£2 = arctan (^ 2 /^ 2 )- 
It can be seen that the phase difference (pi — changes 
over time in a limited range, which means the phase syn¬ 
chronization. Another proof of the existence of phase 
synchronization is a specific form of the phase diagram 
in Fig.4d. In Fig.4e we have shown time-dependence of 
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FIG. 3: Chart of Lyapunov exponents for the system © , a — 
35, c — 25, d = 5, e = 35, k\ — 430, b\ — 12. Parameters & 2 , 
&2 are varied along the Route 1 in Figure 2. Poincare sections 
and phase portraits in the selected points on the Lyapunov 
chart are presented on the lower part of the Figure. 


FIG. 4: Illustrations of synchronous quasiperiodicity, — 
320, /x = 0.7. a), b) - phase portraits of both subsystems, 
c) the time-dependence of the phase difference, d) phase dia¬ 
gram, e) the time-dependence of subsystem amplitudes. 
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the first and second subsystems amplitudes R\ and R 2 . 

The second scenario. Now in the second subsystem we 
will vary the parameters along the Route 2 in Figure 2 . 
The first subsystem generates hyper chaotic oscillations; 
we choose the selected point k\ = 430, 61 = 12 on this 
line. Thus, we are going to discuss the case of interaction 
of two subsystems with hyperchaos. Since every subsys¬ 
tem has two positive exponents in their own Lyapunov 
spectrum, on Lyapunov chart it is logical to expect the 
area HC4 with four positive Lyapunov exponents. This 
idea is confirmed by Fig.5a, that demonstrates regimes 
with two, three and four positive Lyapunov exponents. 
It is remarkable, that at large coupling not only the num¬ 
ber of positive Lyapunov exponents is reduced, but the 
regions of two-frequency quasiperiodicity (although very 
narrow) arise. Thus, the interaction of the hyperchaotic 
subsystems can lead to quasiperiodicity. This statement 
is illustrated in Fig.5b, which presentes an enlarged frag¬ 
ment of the Lyapunov chart. We can see the area of 
quasiperiodic oscillations and built-in system of Arnold 
tongues. 


238 k 2 252 

FIG. 5: The chart of Lyapunov exponents of model © and 
its enlarged fragment. 


Figure 6 is a plot of the largest four Lyapunov expo¬ 
nents along the line = const in Fig.5a. An interesting 
point is a noticeable change in the irregularity of the plots 
in comparison with Fig.l. Thus, hyperchaos arising from 
the interaction of subsystems is rather rough, i.e. weakly 
dependent on parameters. This fact may be important 
for applications. 

Thus, in the present report we deal with the prob¬ 
lems of the interaction of the subsystems that can exhibit 
quasiperiodicity and hyperchaos in autonomous mode. 
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FIG. 6: Plot of the largest Lyapunov exponents as a function 
of the coupling parameter along the line &2 = 285 in Fig.5. 
Phase portraits at g — 0.9 are presented. 


For coupled subsystems with autonomous quasiperiodic¬ 
ity we detect the regime of synchronous quasiperiodicity. 
One of interesting results is the emergence of quasiperi- 
odic oscillations for the interacting hyperchaotic systems. 
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